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Abstract 

We present the algebraic Bethe Ansatz solution for the vertex model recently proposed 
by Zhou H as the classical analog of the Bariev interacting XY chains. The relevant 
commutation rules between the creation fields contain the Hecke symmetry pointed out 
recently by Hikami and Murakami Q. The eigenvalues of the corresponding transfer 
matrix are explicitly given. 
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Some years ago, Bariev has formulated a model of interacting XY chains and solved it 
exactly, in one-dimension, by the coordinate Bethe Ansatz approach Jl|. The model couples 
two XY models and its Hamiltonian on a lattice of length L can be written as [|1], |2|. |5[ 

L 

H = Y, {(vf^+i + °l°t+i) ex P(« r ii^i) + (^\ r +i + n T t+x) ex P(V* + ^ r )} (!) 
1=1 

where af and rf 1 are two commuting sets of Pauli matrices acting on site i, and a is the coupling 
constant. Soon after that, Bariev has generalized this model to include many coupled XY 
chains. After a Jordan- Wigner transformation, the model (1) and its generalizations, can 
be seen as an electronic system where the hopping term depends (asymmetrically) on the 
occupation number of the site itself 0, Q. In this sense, the coordinate Bethe Ansatz solution 
of these models has been used in refs. [|] to explore the finite size behaviour, the excitations 
as well as some related conductivity properties. 

The quantum integrability of (1), however, has only been recently proved by Zhou || in 
terms of the quantum inverse scattering approach. Zhou || was able to construct the two 
dimensional vertex model whose transfer matrix commutes with the Bariev Hamiltonian (1). 
The purpose of this letter is to show that such underlying vertex model can be diagonalized 
by the algebraic Bethe Ansatz |7|, [|. We recall that this method is a powerful mathematical 
technique, which can provide us with information concerning the properties of the Bariev chain 
within a unified perspective. Our formulation is strongly inspired by our recent construction of 
the Bethe states of the Hubbard model by means of the quantum inverse scattering approach 



HI]. We remark, however, that the structure of Zhou's i?-matrix is a bit different than that 



appearing in the Hubbard model |9|, Hj. Indeed, we shall see that the appropriate parametriza- 
tion of the quantum i?-matrix associated to the Bariev chain is different and, in fact, much 
simpler than that we have found for the Hubbard model JTTI . 



The quantum .R-matrix R(X,p) found by Zhou |J has fifteen non-zero Boltzmann weights, 
and following ref. || we denote them by p«(A, p), i = 1, • • • , 15. However, many of the weights 
are related to each other under some functional properties, such as Pi(\, n) = —pj([i, A); Pi(\, p) - 
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Pj(hX, hp); Pi(\, p) = hpj(X, p) []• The parameter h is given in terms of the coupling constant 
a by h = exp(a). For explicit expressions we refer to ref. ||. Here we quote only few extra 
identities which we found relevant in the course of our calculations. These are given by 

pi 5 (X, p)[p 9 (X, p)+p^X, p)] = p 5 (X,p)p G (X,p), Pa{\ p)pi{K P)+Ph{\ p)Pis{\ P) = Pe(X, p)Pi{\ p) 

(2) 

Pi2{\ p)[pv{\ p)+Px{K p)] = p 5 (X,p)p 4 (X,p), p 5 (A, p)pi(X, p)+p 15 (X, p)p 6 (X, p) = p 5 (X,p)p 10 (X,p) 

(3) 

In order to diagonalize the transfer-matrix of the classical vertex model corresponding to 
the Bariev chain we can basically follow the main steps of our recent algebraic construction of 



the Bethe Ansatz for the Hubbard model [11]. We take as the reference state |0), the standard 
ferromagnetic vacuum where all the spins are in the "up" eigenstate of erj and r?. We solve 
the Yang-Baxter algebra for Zhou's .R-matrix by writing the monodromy matrix T(A) in the 
auxiliary space as 

/B(X) B(X) F(X)\ 
T{X)= (7(A) i(A) B*{X) (4) 

\C(X) C*(X) D(X)J 
where -B(A) (B*(X)) and C(X) (C*(X)) are two component vectors with dimensions 1 x 2(2 x 1) 
and 2 x 1(1 x 2), respectively. The operator A(X) is a 2 x 2 matrix and the other remaining 
operators are scalars. The transfer matrix T(A) is the trace of T(X) on the auxiliary space, 
and the eigenvalue problem becomes 

[B(X) + ]T A aa (X) + D(X)] |$„(Ai, • ■ ■ , An)) = A(A, {A,}) |$„(A 1; ■ ■ ■ , A n )) (5) 

a=l 

The set of variables {Ai, • • • , X n } parametrizes the multi-particle Hilbert space by the action 
of the creation fields on the reference state |0). The operators B(X), B*(X) and F(X) play the 
role of creation fields while C(X),C*(X), C(A) and A ab (X), for a ^ b = 1,2, are annihilators. 
This means that the monodromy matrix (4) has a triangular form when acting on the reference 
1 For instance, from ref. [^), it is possible to check that pn(X, /i) = pg(\/h, fi/h); pn(\,ii) — 
P2(\p>)p2(^/h,p/h); p 5 {X,p) = -pi2(p,X). 
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state |0). In addition, we have the following "diagonal" identities 

B(A)|0) = |0), D(\) |0) = [A] 2L |0) , A aa (X)\0) = [Xh} L \0),a=l,2 (6) 

A crucial step in algebraically solving the eigenvalue problem (5) is to find the appropriate 
commutation rules between two fields of B (A) or B* (A) type. Similar to what happens for the 



Hubbard model ||11|| , their commutation rules are equivalent, because they generate as a new 
operator only the common creation field F(X). Remarkably enough, these commutation rules 
already encode the basic underlying hidden symmetry of the Bariev chain || . We can see this, 
for instance, in the commutation relation f\ between the fields B(X) and B(fi) 



B(X) ® B(fi) = B(fi) ® B(X)f(X,fi) + 



-{p 5 (A, n)F(\)B(n) + P5 (fM, X)F(fx)B(X)} (7) 



Pg(A,/i) 

— * 

where the vector £ and the matrix f(A,/x) have the following structures 

(I 0\ 

a(A,/i) 6(A, /i) 

b(X, n) a(A,/i) 

\0 1/ 
and functions a(A,/i), a(X, fi) and &(A,/i) are given by 

A(l-/i 2 ) Ml-^ 2 ) hn , /'(A-/') 



£=(01 /i- 1 0), f(A,/i) 



a(A,/x) 



(9) 



A — h 2 fx X — /i 2 /i ' A — h 2 n 

The structure of the Boltzmann weights of the matrix f(A, /i) are the same of that appearing 
in the 6-vertex model with an azimuthal anisotropy 77 given by 77 = iln(h) = ia. In order to 
see that, we have to introduce the following parametrization 



A = exp[ik(X)} 



(10) 



and consequently the Boltzmann weights can be rewritten in terms of the difference k 
k(X) — k([i) and the anisotropic constant 77 as 



exp(ik/2) sin(r7) 

a ( k ) = — ■ a 10^ \ ' °W 
sm{k/2 + rj) 



exp(— ik/2) sin(i7) 
sin(/c/2 + rj) 



b{k) 



sin(fc/2) 



sin(/c/2 + 77) 



(11) 



We remark that identities such as (2) and (3) are important in the simplification of the commutation rules. 



By means of a transformation which preserves the Yang-Baxter equation, the so-called 
symmetry breaking transformation [13], the Boltzmann weights a(k) and a(k) can be sym- 



metrized in order to give the standard anisotropic 6- vertex model. In fact, the presence of 
the asymmetric version of the 6-vertex model in the commutation rules is a clear sign that 
the underlying symmetry is of Hecke type. It is not difficult to see, from the asymmetric 
vertex model (8), that we can construct a braid operator which appears as the generator of 
the Hecke algebra ||13|| . We remark here that such symmetry was first noticed by Hikami and 
Murakami || in the context of the lattice Schrodinger equation for the "fermionic" formulation 
of the Bariev Hamiltonian (1). The only subtle point is the minus sign on weight b(k). The 
Yang-Baxter equation is invariant under b(k) — > —b(k), and the sign ± can be interpreted as 
periodic/antiperiodic boundary conditions when the size of the quantum Hilbert space is odd 



Now, the construction of the eigenvectors and the eigenvalues goes fairly parallel to the 
formulation we have recently presented for the Hubbard model (Tl| . Basically, we have to adapt 



the " diagonal " commutation rules of ref. (TT| by taking into account the specific weights of 
Zhou's i?-matrix and also consider the convenient parametrization given in equation (10). Here 
we only present our final results for the eigenvalues of the "covering" vertex model. Many other 
results, as well as the main technical steps we have developed will be presented in a separate 
publication JE2] , together with the detailed algebraic solution of the Hubbard model • We 



remark that, recently, the exact expression for the eigenvalue appears to be very important in 
the study of finite temperature properties of integrable models [13, [T3], [TjJ. We found that the 



eigenvalue of the transfer matrix associated to the Bariev chain is given by 

a(a, {a,}) = n '^j^ + a 2l n + ft ^-±^a<» ( a, { A i}) ( i 2) 

i=l z i=l % i=l 1 

where A^(A, {Aj}) is the eigenvalue of the vertex model defined by the auxiliary i?-matrix 

f(X,fi) in the presence of inhomogeneities. Furthermore the variables {Aj} are constrained by 
3 For an even size the sign does not matter. We also recall that in the fermionic formulation of the Bariev 
chain this sign is positive. 
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the Bethe Ansatz equation 

[\ih]~ L = -(-l) n A (1) (A = Xi, {Xj}), i = l,-..,n (13) 

The auxiliary problem can be solved by using the standard 6- vertex formulation of Faddeev 
et al 0, ||, adapted to include the inhomogeneities {Xi}. In the diagonalization procedure, it is 
necessary to introduce the auxiliary variables {pLj} and the eigenvalue A' 1 ) (A, {Xi}, {fAj}) reads 

mi n mi 

a»'(a, {a.}, w) = n w + n ka, a.) n ra («> 

where the variables {fij} satisfy the equation 

fi%,>A,) = -n^4, j = l,---,m (15) 

t=l k=l b ^k,fij) 

Finally, all these results can be combined in order to give us the eigenvalue and the Bethe 
Ansatz equations. At this point, to cast the final results in a convenient form, we redefine the 
variables A, {Aj}, and {fij} by 

Aj/i = exp(iki) , fij = exp(iAj), A = exp(ifc) (16) 

In terms of these new parameters the expression for the eigenvalue is 

A(*. {h}, {A,}) = ft ™W 2 + W-'>M + exp (i2Lk) ft + 

' 1 X ' 1 J/; z sin(^/2 - fc/2 + 77/2) P 1 ; t=x 1 - h/2 + rj/2) 

+ exo U(k - n)L] I ft i™W* + W-vM ft sin(A J /2 - fc/2 + g) 
+ exp r^Lj |U s . n( ^ /2 _ k/2 + ^ /2) n s . n(A ^ /2 _ k/2) + 

" z cos(A:/2 + h/2 - r//2) ™ sin(A;/2 - A./2 + g) | 

sin(fc/2 - A; 4 /2 + 77/2) /J sin(fc/2 - A.-/2) J 1 j 

and the nested Bethe Ansatz equations are given by 

exv(ik-L) - -(-1)™ ft ^(^2-^/2 + ^2) _ 

exp(^L)- ( 1) J| sin(fei/2 — Aj/2 — rj/2) ' ' " (18) 

" sin(A J /2 - k/2 - r//2) _ " sin(A J /2 - A fc /2 - 77) 
1 J sin(A j /2 -fci/2 + 77/2) " ^ sin(A,/2 - A fc /2 + 7/) ' J ' ' 1 ' 
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This last equation is similar to that found early by Bariev [p], @] in the context of the 
coordinate Bethe Ansatz approach, as should be. In order to recover the results of Bariev 
[[I], for the eigenenergies E(L) of the Hamiltonian (1), we just have to take the logarithmic 
derivative of the transfer matrix eigenvalue at the point A = 0. This calculation leads us to 

L 

E{L) = 2/i^cosOi) (20) 
i=i 

We conclude this letter with the following remarks. The extra signs we have found in 
the Bethe Ansatz equations (18,19) are typical of the " bosonic " formulation (1) of the 
Bariev chain. They can be related with peculiar boundary conditions [JTTl [H]], and they are 
not present if one formulates the diagonalization problem for the "fermionic " version of (1). 
Our algebraic formulation has an invariance under h — > h , which is in accordance with 
the symmetry of the Bariev chain ( a — > —a) [1X1] . Following the results of Shiroischi and 



Wadati ||18| , there exists a way of generating a generalized Bariev chain from Zhou's R- matrix. 
Defining the vertex operator |18| £ e °(X) = PR(X,9q), where P is permutator and R(X,9q) is 
Zhou's i?-matrix, we can define a one-parameter (#0) family of vertex models by the transfer 
matrix T e °(A) = Tr a [C e ° L {\) ■ ■ ■ C e a \{\)). Such vertex model can be diagonalized following 
the basic steps we presented so far. The main change is concerned with the action of the " 
diagonal " operators on the reference state |0). In this case we find B(X)\0) = [pi(A, 6* )] L |0), 
A aa (X)\0) = [p 3 (X,6 )} L \0), D(X)\0) = [p 9 (AA)] L |0), in such way that only the terms which 
are proportional to the power of L change in the expressions (12) and (13). Finally, since 
the parametrization (10) is quite simple, it seems interesting to re-investigate the Yangian 
symmetry as well as the analytical properties of the transfer matrix associated to the Bariev 



chain in light of the recent results of refs. [19, 20 
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